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Abstract 



The extremely small probability of quantum tunneling through an almost 
classical potential barrier may become not small under the action of the spe- 
cially adapted nonstationary field. The tunneling rate has a sharp peak as a 
function of the particle energy when it is close to the certain resonant value 
defined by the nonstationary field (Euclidean resonance) . Alpha decay of nu- 
clei has a small probability since the alpha particle should tunnel through a 
very nontransparent nuclear Coulomb barrier. The incident proton, due to 
the Coulomb interaction with the tunneling alpha particle, plays the role of 
a nonstationary field which may result in Euclidean resonance in tunneling 
of the alpha particle. At the resonant proton energy, which is of the order of 
0.2 MeV, the alpha particle escapes the nucleus and goes to infinity with no 
influence of the nuclear Coulomb barrier. The process is inelastic since the 
alpha particle releases energy and the proton gains it. This stimulation of 
alpha decay by a proton constitutes a new type of nuclear reaction. 



PACS number(s): 03.65.Sq, 42.50.Hz, 24.10.-i 
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I. INTRODUCTION 



A control of processes of quantum tunneling through potential barriers by external sig- 
nals is a part of the field called quantum control which is actively developed now, see, for 
example, Ref. [1] and references therein. Excitation of molecules, when one should excite 
only particular chemical bonds [2-4] , formation of programmable atomic wave packets [5] , a 
control of electron states in heterostructures [6], and a control of photocurrent in semicon- 
ductors [7], are typical examples of control by laser pulses. A control of quantum tunneling 
through potential barriers is also a matter of interest, since tunneling is a part of many 
processes in nature. The computation of probability for a classically forbidden region has 
a certain peculiarity from the mathematical stand point: there necessarily arises here the 
concept of motion in imaginary time or along a complex trajectory [8-10]. The famous semi- 
classical approach of Wentzel, Kramers, and Brillouin (WKB) [8] for tunneling probability 
can be easily reformulated in terms of classical trajectories in complex time. The method of 
complex trajectories is also applicable to a nonstationary case [11,12]. The method has been 
further developed in papers [13-17], where singularities of the trajectories in the complex 
plane were accounted for an arbitrary potential barrier (see also [18]). Recent achievements 
in the semiclassical theory are presented in Refs. [19-23]. 

Let us focus on the main aspects of tunneling under nonstationary conditions. When the 
electric field £ cos VLt acts on a tunneling particle of the initial energy E, it can absorb the 
quantum Ml (with the probability proportional to the small parameter S 2 ) and tunnel after 
that in a more transparent part of the barrier with the higher energy E + Ml. The pay in 
the absorption probability may be compensated by the probability gain in tunneling. In this 
case the system tends to absorb further quanta to increase the total probability of passing 
the barrier. This mechanism of barrier penetration is called photon-assisted tunneling. If 
Ml is not big, the process of tunneling, with the simultaneous multiquanta absorption, can 
be described in a semiclassical way by the method of classical trajectories in the complex 
time [13-17]. When a tunneling particle of the energy E is acted by a short-time pulse, the 
tunneling probability is associated with the particle density carrying away in the outgoing 
wave packet. The particle energy after escape is E + 5E, where the energy gain 8E = Ntrw, 
should be extremized with respect to the number of absorbed quanta N and the energy hio 
of each quantum [24] . 

According to the perturbation theory, the both probabilities, of absorption and emission 
of quantum, are small being proportional to S 2 . After absorption, SE is positive which 
enhances the total probability due to increase of the tunneling rate. After emission, SE is 
negative and the particle should tunnel with a smaller energy in a less transparent part of 
the barrier; in this case there is no gain in the probability due to tunneling as for absorption. 
At first sight, tunneling in a nonstationary field cannot be assisted by an emission of quanta 
of this field since one should pay in probability twice. Thus, one can expect the double loss 
in probability due to the emission of quanta and the reduction of tunneling transparency. 
This conclusion is correct as soon as the nonstationary field is small and the perturbation 
theory is applicable. 

Under increase of the nonstationary field the process of tunneling with emission of quanta 
becomes completely different compared to one, predicted on the basis of the perturbation 
theory [25]. The crucial role here plays the fact, that the non-perturbative wave function is 
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mainly determined by its big phase (real or imaginary). In this case there is no the double 
loss in probability. The quantum process is not simply reduced to separate emission and 
tunneling. However, the phase behavior can be interpreted in the way of enhancement of the 
total probability due to emission processes which competes now with the reduction of the 
tunneling transparency. The competition between the enhancement of the total probability 
due to emission and the reduction of it due to tunneling results in the unexpected effect: the 
total probability (defined as a particle density, carrying away by the outgoing wave packet) 
becomes not exponentially small for the certain particle energy Er in the well. This energy 
depends on parameters of the nonstationary field and in a vicinity of E R the probability 
sharply peaks as a function of energy. This reminds, formally, a resonant behavior and is 
called Euclidean resonance [25]. The energy Er can be called the resonant energy. In III 
this phenomenon is deduced from the analysis of the quantum mechanical phases using only 
simple physical arguments. 

As well known, nuclear Coulomb barriers may be very significant in nuclear physics, 
playing a role of blockade for particle approach or escape [26,27]. The famous example of 
such a nuclear process is alpha decay of nuclei which has a small probability since the alpha 
particle should penetrate through a very non-transparent nuclear Coulomb barrier [26,27]. 
As any tunneling process, alpha decay can be influenced by a nonstationary field. The 
duration of artificially generated pulses (see, for example, [28]) is too long compared to the 
nucleus characteristic time of 1CT 21 s which makes impossible their influence on alpha decay. 
A role of a nonstationary field can be played by a moving charged particle which collides the 
decaying nucleus. For example, a proton of the energy of 1 Mev sweeps the tunneling region 
in, approximately, 10~ 21 s. An incident proton, moving towards a tunneling alpha particle, 
reduces its energy due to the Coulomb interaction between them. Since the energy is lost 
but not gained during tunneling, Euclidean resonance may be expected. As shown below, 
this happens. At the certain energy of the incident proton the simultaneous tunneling of 
alpha particle has not an exponentially small probability, according to Euclidean resonance. 

In other words, due to interaction with an incident proton of the certain energy, the 
alpha particle escapes the nucleus and goes to infinity with no influence of the Coulomb 
barrier. This stimulation of alpha decay by a proton constitutes the new type of nuclear 
reaction. Normally, nuclear reactions involve strong interaction at nuclear distances. In 
this case, only the short range start of the alpha particle from the nucleus is due to strong 
interaction and the main physics occurs further, at the larger distance where strong forces do 
not act and only Coulomb effects are involved. This new type of nuclear reactions has a very 
resonant character with respect to an energy of the incident proton. As shown below, the 
typical resonance energy of a proton is in the range of 0.2 MeV. One should emphasize, that 
this resonance results from solely Coulomb effects in contrast, for example, to resonances in 
nuclear physics due to formation of compound nuclei. 

II. PHOTON-ASSISTED TUNNELING 

A penetration of a particle through a potential barrier is forbidden in classical mechanics. 
Only due to quantum effects the probability of passing across a barrier becomes finite and 
it can be calculated on the basis of WKB approach, which is also called the semiclassical 
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theory. The transition probability through the barrier, shown in Fig. 1, is 

W ~exp[-A (E)} (1) 

where 

ME) = \J dx^2m [V(x) - E] (2) 

is the classical action measured in units of h. The integration goes under the barrier between 
two classical turning points where V(x) = E. One can use the general estimate A ~ V/hw, 
where V is the barrier height and u is the frequency of classical oscillations in the potential 
well. A semiclassical barrier relates to a big value V/ftiv ^> 1 . 

What happens when the static potential barrier V(x) is acted by a weak nonstationary 
electric field £(t)7 In this case there are two possibilities for barrier penetration: (i) the 
conventional tunneling, which is not affected by £(t), shown by the dashed line in Fig. 2(a), 
and (ii) an absorption of the quantum Ml of the field £{t) and subsequent tunneling with 
the new energy E + Ml. The latter process is called photon-assisted tunneling. The total 
probability of penetration across the barrier can be schematically written as a sum of two 
probabilities 

V\ fa£a\ 2 ( V -Ml 



where £n is the Fourier component of the field £(t). The second term in Eq. (3) relates to 
photon-assisted tunneling and it is a product of probabilities of two quantum mechanical 
processes: absorption of the quantum Ml and tunneling through the reduced barrier V — Ml. 
The length a is a typical barrier extension in space. When the frequency is high Q > u, 
the second term starts to dominate at sufficiently small nonstationary field (aSn/h) 2 > 
exp(— Q/lj). This is a feature of tunneling processes. Normally a nonstationary field starts 
to dominate at bigger amplitudes (aS^/h) 2 > 1. When the second term in Eq. (3) exceeds 
the first one, further orders of perturbation theory should be accounted which correspond 
to the multiple absorption, shown in Fig. 2(b). 

Let us specify a shape of a field pulse in the form 

m = w 

with the Fourier component £q = n£9 exp(— 9\Q\). In this case, in addition to the steady 
particle flux from the barrier, an outgoing wave packet is created which carries away the 
certain particle density. Then the probability W of the transition through the process of 
absorption of iV quanta and subsequent tunneling with the higher energy E + NMl, shown 
in Fig. 2(b), is 

W ~ (^Y^j exp[-A (E + Nhfl)] = (^9a£y N exp(-A) (5) 



where 



29 

A = — 5E + A (E + 5E) (6) 
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Here the total energy transfer SE = NhQ is introduced (ft > 0). The maximum transition 
probability through the barrier is determined by some finite value of SE, which provides 
a minimum of A and is defined by the condition dA(E + SE)/dSE = 0. An existence of 
such a minimum is possible if a grow up of small SE reduces A, that is, under the condition 
29 < h | dA (E)dE | of sufficiently short pulses. In other words, sufficiently short and not 
very small pulses (however, still much smaller than the static barrier field) strongly enhance 
tunneling by photon assistance. 

Eq. (5) omits some details and it is rather an illustration of a tunneling mechanism with 
quanta absorption. For example, the accurate perturbation theory starts with the linear 
£-term. The exact nonperturbative theory, which is a generalization of the conventional 
semiclassical approach, results in the same A. The semiclassical approach is also sensitive 
to the sign of £ since the classical energy transfer is determined by £{t) dx/dt, where x(t) is 
a classical particle trajectory. The strong photon-assisted tunneling exists only at positive 
£. 



III. EUCLIDEAN RESONANCE 



Besides the absorption of quanta in Fig. 2, also the emission is possible, shown in Fig. 3. 
This process is provided by negative frequencies Q < and the energy transfer is also 
negative SE < 0. In this case Eq. (5) gives A = 20\SE\/h + A (E — \SE\) which does not 
correspond to any extreme since dA/d\SE\ > 0. So, on the basis of perturbation theory, 
one can conclude that emission processes cannot enhance a barrier penetration at least for 
the pulse amplitude on the border of applicability of the perturbation theory a£ < h/6. 
May an emission process enhance a tunneling rate when the signal amplitude is not small 
h/6 <a£? 



A. Phase connection 

At the big pulse amplitude £ its contribution to the classical action, generally speaking, 
is big, compared to Planck's constant, 

a£d > h (7) 

and the wave function ip ~ exp(ix) is mainly determined by its big phase x, which can be 
imaginary as well. Another condition of a big phase is a slow varying pulse 

£ « 6 (8) 

since a big phase should be built up during a long time. Here V is the barrier height and, 
therefore, h/V is some intrinsic time of the problem. One has to note, that within the 
conditions (7) and (8) the pulse amplitude can be still less compared to the static field of 
the barrier V/a. 

Suppose %i and x 2 to relate to the classical turning points (1) and (2) in Fig. 3. At 
t — > ±oo, when £(t) = 0, there is a conventional tunneling through the barrier. For 
a symmetric in time pulse £{t) the modulus of the wave function of the outgoing particle 
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\ijj(x2, t) | has an extreme (maximum) value at t — which relates to the maximum amplitude 
of the outgoing wave packet at t > 0. According to Feynman [29], an extreme wave function 
corresponds to a classical trajectory of the particle connecting the points {xi, 0} and {x 2 , 0}. 
But in the present situation there are no classical trajectories under a barrier. 

However, one can find a connection between two constants ip(xi, 0) and ip(x2, 0) without 
an exact solution of Schrodinger equation. The main point of this procedure is the possibility 
to define the wave function mainly by its big phase (real or imaginary) which is true under 
the conditions (7) and (8). Then one can consider a particle state at any coordinate as 
in classical mechanics. Since the moment t — corresponds to the extreme situation of 
maximum output, one can look for an alternative extreme way to connect phases of ip(x±, 0) 
and ip(x2, 0). Let us find some formal path from (1) to (2) which relates to an extreme total 
phase with respect to the energy of the final state (2). 

Suppose, that the particle state (i) at x = Xi is chosen with the same energy E — \SE\ 
as at the point (2) in Fig. 4. The particle from the position (i) can tunnel to the position 
(2) and the acquired phase is imaginary 



^(x 2 ,0) ~ ij(xi,0)exp 



--A (E-\5E\) 



(9) 



On the other hand, the particle can absorb N - 
to the state (1) where tp(xi,0) ~ -^(^(J) (£ n ) N . 
phase 



: \SE\/hfl quanta, as in Fig. 4(a), and go 
This also leads to an imaginary acquired 



ip(x!, 0) ~ ipfa, 0) exp f -^1^1 



(10) 



We use here the expression for Sq. Analogously, the particle can emit iV quanta and go from 
(1) to (i), as in Fig. 4(b). This leads to another connection of (1) and (i) 



i/j(xi, 0) ~ t/>(xi, 0) exp ( 



(11) 



According to Eqs. (9) - (11), the connection between (1) and (2) can be written in the form 



ip(x 2: 0) ~ exp 



A Ao (E-\5E\ 



c„ exp ( ^ ()/•; J + r,,e.\p I 



^(xi,0) 



;i2) 



where c a and Cj, are constants. Only the first term in Eq. (12) provides an extreme of the 
total phase and, therefore, one should put c a — 1 and q, = 0. This is an alternative extreme 
way of phase connection. The relation (12) corresponds to a formal (no direct analogy with 
the exact solution) path connecting the phases (1) and (2). The extreme (at the moment 
t = 0) transition probability from (1) to (2) is determined by the found extreme phase 
difference in (12) 



W(l -> 2) 



^2,0) 

ip(xi,0) 



exp [-A(\SE\)\ 



(13) 
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where 



A(\SE\) = A (E- \6E\) - — \5E\ (14) 

As also follows from the solution of Schrodinger equation, weakly depends on the 

pulse amplitude S under the conditions (7) and (8). The optimum energy transfer \8E \ 
should be determined from the extreme condition dA(\SE\) / d\SE\ = 0, which reads 

29 dA (E-\5E\) 

T = dE (15) 

and determines the extreme value A = A(\5E \). Again, the physical transition probability 
at t = 0, W ~ exp(— A), can be considered either as an extreme in t of the exact solution 
or as an extreme in \SE\ in the above phase connection procedure. 

The phase connection, illustrated in Fig. 4 is applicable for big pulses K/6 < a£ when the 
wave function is mainly defined by its big phase (real or imaginary) and the phase difference 
between (1) and (i) is of the opposite sign compared to one between (i) and (1). This phase 
connection, as follows from the solution of Schrodinger equation, is applicable in our case 
of decay of a metastable state. But it does not mean, that its applicability is automatically 
valid in other situations with a big phase. For example, a penetration of an incident particle 
into the potential well is not described by this method. 



B. Resonance conditions 



An applicability of the phase connection, besides the conditions (7) and (8), is restricted 
by the inequality 

exp(-A) < 1 (16) 

A pulse amplitude S should be negative, since the classical energy transfer is determined 
through the classical trajectory x(t) as S(t)dx/dt. A is determined by the particle energy 
E in the well and by the pulse duration 9. Suppose 9 to be fixed and the particle energy E 
to vary. Then the optimum energy transfer \SE (E)\, determined by Eq. (15), is a function 
of E and one can define the certain energy E R {9) as a solution of the equation 

29 

A [E - \6E (E)\] - — \5E (E)\ = (17) 

At the particle energy close to E R) A ~ 29[E R — E]/h and the peak in tunneling probability 
at t — 0, related to the particle density in the outgoing wave packet, is 

W~exp(-™[E R (0)-E]) (18) 

The formal applicability of these relations is W 1, nevertheless, the probability peak 
can dramatically grow up (upon approaching E R ), for example, from 1CT 37 to 1CT 2 . The 
energy E R has the order of magnitude of the barrier height V and, hence, 29(E R — E)/h ~ 
(9V/h)(E R — E)/E R . This means, that the tunneling probability at t — has a sharp peak, 
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like a resonance, as a function of the particle energy E near Er. The effect may be called 
Euclidean resonance when Er plays a role of the resonant energy. The origin of the word 
"Euclidean" is explained below. 

Euclidean resonance also can be treated in another way. Suppose the energy level in the 
well to be fixed. Then one can adjust the pulse parameter 9 to meet the condition E R = E 
when the barrier becomes almost transparent at t — for the energy E. Note, that for other 
particle energies the barrier remains low-transparent. 

In contrast to photon-assisted tunneling, which has a connection with the perturbation 
theory (see Fig. 2), the phenomenon of Euclidean resonance is completely non-perturbative. 
One can construct some analogy of Euclidean resonance by adding the narrow potential well 
— v(b)5(x — b) to the potential V(x) in Fig. 1. The distance b is within the under-barrier 
region and the positive coefficient v(b) is chosen to get the same energy level E in the 5- well 
as in the main well (resonance tunneling). Suppose b to be close to the main well in Fig. 1 
and the peak of the wave function at x = b is of the same order as in the main well. Let 
us move b slowly (compared to the time h/V) away from the main well to get the big peak 
of the wave function at some point under the barrier. After this v(b) is switched off fast 
and the remaining distribution, which is not exponentially small, goes partly outside the 
barrier. So, this non-stationary mechanism provides the outgoing wave packet which is not 
exponentially small. As in Euclidean resonance, two issues are important in this example: 
(i) the nonstationary potential should be not small and (ii) it should be very precisely chosen 
to get the condition of resonance tunneling, otherwise, the peak of the wave function inside 
the barrier would be exponentially small. 



Let us consider an electron emission from a metal, left in Fig. 5, to the vacuum due to 
the applied electric field £ + £(t) where £ is a constant. The energy E is supposed not 
to be above the Fermi level. An electron emission occurs by tunneling through the barrier 
V(x) — x£(t) where the pulse acts only at x > 0. The static barrier is V(x) = V — x£ at 
x > and V(x) = at x < 0. 

The conventional WKB action (2) has the form 



IV. AN EXAMPLE OF EUCLIDEAN RESONANCE 



A (E) 



^(V-E)roo(E) 



(19) 



where 



roo(E) 



T y/2m(V-E) 



(20) 



The relation (15) turns to t 00 (E — \SE\) 



- 9 and the optimum energy transfer is 



5E, 




V + E 



(21) 



o 



2m 



Eq. (17) defines the resonance energy 



E R {9) = V 



9^£l 



(22) 



6m 
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One can note, that, under the resonance condition, the input energy E R is connected to the 
output one E R - \8E\ as 3(V - E R ) — V— (E R - \5E\). 

These formulae relate to the Lorentzian shape of a pulse (4). What happens for the 
nonstationary field £{t) = £ cos fit or of some other shape? The above arguments, based on 
general physical principles, are not sufficient to answer this question and a more sophisticated 
treatment is required. This is considered in V. 



V. TRAJECTORIES IN IMAGINARY TIME 

According to Feynman [29] , when the phase of a wave function is big, it can be expressed 
through classical trajectories of the particle. But in our case there are no conventional 
trajectories since a classical motion is forbidden under a barrier. Suppose a classical particle 
to move in the region to the right of the classical turning point (2) in Fig. 5 and to reach 
the point (2) at t — 0. Then, close to the point (2), x(t) = x 2 + ct 2 (c > 0) and there is no a 
barrier penetration as at all times x(t) > x 2 - Nevertheless, if t is formally imaginary, t = it, 
the penetration becomes possible since x(ir) = x 2 — cr 2 is less then x 2 . Therefore, one can 
use classical trajectories in imaginary time to apply Feynman's method to tunneling. 



A. Newton's equation and classical action 

A classical trajectory has to satisfy Newton's equation in imaginary time 

d 2 x dV(x) 

where V(x) is the static barrier in Fig. 5. The classical turning point (2) in Fig. 5 is reached 
at t = with the initial condition 



dx 



point (2) (24) 



T=0 



The point (1) is reached at t — ir when 

x (iT ) = point (1) (25) 

The "time" To is expressed through the particle energy E before the barrier where the 
nonstationary field does not act 

The conditions (24), (25), and (26) define the solution x(ir) of Eq. (23). This solution, in 
turn, defines the extreme (in time) transition probability related to the particle density in 
the outgoing wave packet 

W ~ exp(-A) (27) 
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where A is the classical action in units of Planck's constant 



TO 



dr 



y(^) +V(x)-x£{ir) 



E 



(28) 



r can be treated as "time" of motion under a barrier. Eq. (28) holds for any barrier which 
is zero at x < 0. Without a nonstationary pulse, £ = 0, tq coincides with its static value 



TOO - \lj J g_ g 



(29) 



and the action (28) turns to A (E) (2). The integration in Eq. (29) goes between classical 
turning points where V(x) = E. According to classical mechanics, 



2r, 



oo 



dA (E) 
dE 



(30) 



B. An example of imaginary trajectories 

Let us consider the particular pulse (4). In imaginary time S(ir) = 1/(1 — t 2 /8 2 ) diverges 
at r = 9 and this sets r in Eq. (28) close to 9. A difference between r and 9 is small at 
small pulse amplitude a£ <C V. For this reason, the "time" interval (t — 9) near r , when 
S(ir) is not small, weakly contributes to the action. During the short "time" (r — 9) the 
particle energy reduces (because £ < 0) down to (E — \SE\), so that 



too(E - \6E\) ~ 6 
and the particle motion at < r < 9 can be considered to be free 



h 



m ( dx 



+ V(x) - E 



(31) 



(32) 



Adding and subtracting the energy transfer \5E\ in the right-hand side of Eq. (32) and using 
the equation (31), one can arrive at 



A = Aq{E-\SE\ 



29 
1 



5E\ 



(33) 



With account of Eqs. (30) and (31), it is obvious, that Eqs. (14) and (33) give the same 
result. 



C. Why a small pulse may enhance tunneling 

It is already been shown in a simple way in II, that the effective nonstationary field 
in tunneling processes enhances compared to other ones. This enhancement is clearly seen 
in the method of trajectories in imaginary time where an influence of a pulse becomes 
substantial under the approximate condition £{ir) ~ £q. Here £ ~ V/a is the field of a 
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static barrier. In the case of the pulse (4), £{ir) ~ £Q/{j — 0) is enhanced compared to £ 
since one can choose r very close to 6. But this enhancement is not unlimited due to the 
semiclassical condition of slow varying in time (r — 6) <C h/V. Thus, the condition on the 
pulse amplitude coincides with (7) or it can be written as 

t < T (34) 

Since Aq ~ V9/H is big, tunneling can be strongly influenced by a pulse amplitude £ which is 
less then the static field of the barrier. This statement is true for more general forms of non- 
stationary fields. For example, the monochromatic £ cos fit and the Gaussian £exp(— fl 2 t 2 ) 
pulses become exponentially big in imaginary time and their stimulation of tunneling occurs 
at small amplitudes £ < £ as well. 



D. Euclidean resonance in real time 

The metric in relativity x 2 + y 2 + z 2 — c 2 t 2 is Euclidean in imaginary time t = ir, since 
the all coefficient become positive, and the action (28) or (32) is called Euclidean action. 
Extending this analogy, one can name the present phenomenon Euclidean resonance. 

The above approach is valid when exp(— A) <C 1, but upon reduction of A this condition 
may be violated. When exp(— A) becomes no small one should account further contributions 
exp(— 2A), exp(— 3A), etc. This is equivalent to an account of multi-instanton contributions. 

For the particular barrier V(x) = V — x£ in presence of a non- stationary field one can 
build up a bridge between trajectories in imaginary time and the solution of Schrddinger 
equation in real time. One can find this solution in the form ip(x, t) = a(x.t) exp[iS(x, t)/h], 
where S(x, t) is the classical action and one can obtain any correction in h in the pre-exponent 
a(x,t). The result for the imaginary part of the action is shown schematically in Fig. 6. 
Without a pulse, the solution to the left of the point "exit" , vp = C\ipi + c^i consists of 
the dominant, ipi ~ exp(iSi/h), and the sub-dominant, ip 2 ~ exp(iS 2 /h), solutions. ip 2 ( x = 
0) ~ tpi(x — 0) exp(— Aq) is exponentially small, at the point "exit" ipi ~ ip2 ~ exp(— A Q /2), 
and to the right of "exit" there is only an outgoing wave of the amplitude exp(— A /2). This 
is a picture of decay of the metastable state, localized near x — 0, through a static barrier. 

When £(t) is not zero, the third solution, ip 3 ~ exp(iS 3 /H), appears which is shown by 
the curve (3) in Fig. 6. The maximum of this solution 

if> 3 Mt),t]~c 3 exp(-A/2) (35) 

(A is defined by Eq. (28)) is reached at the classical trajectory x d (t) in real time, when 
x c ;(±oo) = oo and x c ;(0) is the minimum value. We omit not strong effects of quantum 
smearing of the wave packet. At a fixed moment of time t < one can find ip 3 (x,t) for all 
x using the generalized semiclassical approach which nowhere breaks down upon sweeping 
over all x. In this case, ip 3 is an independent solution and one should put c 3 = at t < 
since there is no incoming wave. Close to the moment t — the solution tp 2 and ip 3 get 
a tendency to merge at some point, circled by the dashed curve in Fig. 6. In this region 
d 2 S/dx 2 , calculated semiclassically, becomes big and the semiclassical approximation breaks 
down in a vicinity of the circled point. This means, that within the short (non-semiclassical) 
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time interval t ~ H/V the solution ^3 is formed. Then, at t > 0, the semiclassical solution 
recovers at all a; again, but now C3 is not zero, which relates to an outgoing particle. From 
a semiclassical point of view, there is a jump (since a semiclassical approximation does not 
resolve a short time) from zero to one of the coefficient at the third solution at the moment 
t = 0. 

Despite of that the above solution is obtained analytically for the barrier V{x) = V — x£ , 
the presented scenario of stimulation of tunneling by a nonstationary field holds qualitatively 
for a general semiclassical barrier. The jump of the coefficient occurs both in photon-assisted 
tunneling and in Euclidean resonance. 



VI. EUCLIDEAN RESONANCE AND NUCLEAR REACTIONS 

In VI we apply the developed ideas of stimulation of tunneling to nuclear reactions 
where tunneling through a Coulomb barrier is a substantial part of the process. The role of 
nonstationary field is played now by a charged incident particle. 



A. Alpha decay of nuclei 

According to Gamov [26,27], alpha decay of nuclei, is described by tunneling of alpha 
particle through the Coulomb barrier. The potential energy, as a function of the distance R 
between alpha particle and the nucleus, is shown in Fig. 7, where the Coulomb tail olm/R 
(«m = 2(^o — 2)e 2 ) sharply drops at the nuclear size xq. For the alpha decay 

1 5 U - ^Th + a (36) 

z = 92, E = 4.678 Mev, and, according to the liquid drop model, 

x Q = 1.2 [(231) 1/3 + 4 1/3 ] x l(T 13 cm ~ 0.92 x l(T 12 cm (37) 

The WKB tunneling rate is 

W ~ exp [—A M (E, L M )\ (38) 
where the action in units of Planck's constant is 



V8M [ R * Jn /a M , L 2 M 



= Vi» l¥ + 2S- J (39) 

M is the mass of alpha particle and the classical exit point R e is determined by zero of the 
square root. Lm h is the angular momentum. Further we consider parameters under the 
condition 



L 2 M olm 



"0 

The action can be written in the form 



2Mxl x Q 



> E (40) 



, . na M /2M 



K V OL M 



(41) 
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where 



f(p) = VTT^-^\n\i + 2 P + 2VrtT+p)] ; v = l2m,2Mx1 (42) 

p is the ratio of the centrifugal energy and the Coulomb one at the nucleus radius. The 
imaginary time of motion under the barrier 



1 dA M (E } L M ) na M M 

™ (E) = -2 — m — "2^V2E (43) 

weakly depends on the angular momentum Lm under conditions (40). 

According to the semiclassical applicability, the correction to Am due to the second term 
in Eq. (41) should be much bigger then one, that is 



h 2 

< v^o (44) 



32Ma M 



This coincides with the conventional WKB condition in Coulomb field. In this problem of 
alpha decay the condition (44) reads as 0.01 1. Since the parameter p ~ 1, the angular 
momentum is big 



/^>1 (45) 



h V h 2 x 

As one can see, semiclassical conditions are fulfilled well for alpha decay. 



B. An incident proton 

What happens to alpha decay when a charged particle (proton, for example) is stopped 
by the Coulomb field of the nucleus which is ready to emit alpha particle? The classical 
motion of the proton in the Coulomb field of uranium nucleus is described by the equation 

ml Lt, a • 



drle- -P^ - — (46) 
2 J V 2mr r J 

where m is proton mass, a m = zoe 2 , L m is the proton angular momentum, and e is the 
proton energy. The classical trajectory is shown in Fig. 8, where the shortest distance r e 
between the proton and the nucleus (the classical turning point) is given by the zero of the 
square root in Eq. (46). At r < r e the time t in Eq. (46) becomes imaginary and, starting 
at the point r e , the proton reaches the nucleus (the dashed line in Fig. 8) at the "moment" 
t = iT m . Analogously to Eq. (43), the expression for r m is 

na m fm 

Tm{£) = 2fe V27 (47) 
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C. Alpha particle meets proton 



When the uranium nucleus emits the alpha particle, the additional interaction energy 



Vi 



®-int 



int 



\R(ir) 



(48) 



r[iT) 



where a int = 2e 2 , results in a connection of motions of the alpha particle and the proton. 
Now there is a cooperative motion of two particles in imaginary time which starts at r = 0, 
with zero radial velocities, and terminates at the nucleus at the certain under-barrier time To- 
The total energy (E+e) of two particles conserves, e is the energy of the incident proton and 
E is the energy of the alpha particle close to the nucleus at r = r . The interaction energy is 
always small excepting a narrow vicinity of the moment r = t when two particles are close 
to the nucleus where \R — f| ~ xq. In the small vicinity of To the interaction redistributes 
energies, so that the alpha particle leaves the interaction region with the energy E — \SE\ 
and the proton energy becomes e + \SE\. The major part of the interval (to, 0), excepting a 
small vicinity of r , the particles move independently with the redistributed energies, reach 
the point r = 0, and go to infinity in real time having the same energies E — \5E\ and 
e + \SE\. The interaction (48) contributes weakly to the action since it works during a short 
time. 

The both trajectories are shown in Fig. 9, where the exit points are R e = ckm/(E — \ SE\) 
and r e = a m (e + \SE\). As x is small, the both curves are about to merge at r = r , 
otherwise the interaction at this region is not effective. This requires the condition 



T M (E-\6E\) = T m (e + \6E\)~T 
With Eqs. (43) and (47) the condition (49) reads 



(49) 



e+\5E\ 
E — \5E\ 



Ma 2 M 



1/3 



(50) 



D. Cooperative motion of alpha particle and proton 

The cooperative motion of the alpha particle and the proton relates to the Euclidean 
action 



2 f T " 
A = - / dr 




O-M m ( 0f\ Orn aint 

R 2 [dr] r \R-f\ 



E-e 



(51) 



The classical trajectory satisfies Newton's equation resulting from a minimization of A with 
the initial conditions 



dR r 



dr 



dr x 
dr 



OR 

; Ry(0) = r„(0) = ; MR X 



dr y 

L M ; mr x — 



— L m 
(52) 
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where Lm and L m are the angular momenta of particles which conserve for all r excepting 
a close vicinity of r . The another condition is 

R ( T o) = r ( r o) = x (53) 
and t relates to the energy of alpha particle 



E = 



M ( dR\ 2 a M 



(54) 



J T 



In Eqs. (51)-(54) the vectors are defined as R = {R x ,iR y } and f = {r x ,ir y } since in 
imaginary time y-components are imaginary. This can be seen in the proton motion in the 
Coulomb field of the nucleus (with no alpha particle) under the condition analogous to (40) 
for proton 



ma m J \ ma m J ' ' V 2me e ' \ 2eL 2 m/ 

The dependence r(ir) is given by Eq. (46) with t = it. As follows from Eq. (55), r 2 —r 2 = r 2 . 
Fall of the proton to the nucleus means that r(ir) = 0, but r x (ir) and r x (ir) separately are 
not zero. 

In contrast to the big energy transfer \5E\, the interaction (48) results in a very small 
transfer of angular momentum 5L m between alpha particle and proton. On the basis of the 
classical relation 

dS L m a int 



Or \R-r\t 



{r x R y - r y R x ) (56) 



with the estimate R x ~ r x ~ \R — f| ~ a m /e and Eq. (55) for the y-component, one can 
easily deduce, that 5L m /L m ~ ai nt /am 1. As shown below, the energy transfer \5E\ is 
not small since it is determined by by small \R — r\ ~ Xq. The proton can change its angular 
momentum solely by interaction with the nucleus. We consider only spherically symmetric 
nuclei and the proton angular momentum L m conserves. The angular momentum Lm of the 
alpha particle is determined by the interaction with the nucleus and conserves for all r. 

Now it is possible to calculate A in Eq. (51). A contribution to A from the interaction 
part Vi n t is small, almost all "time" r the alpha particle conserves its energy E — \SE\, and 
the same is true for the proton with the energy e — \SE\. For this reason, one can write A 
in Eq. (51) as a sum of two actions of free particles 

A = A M (E - \5E\,L M ) + A m (e + \5E\,L m ) (57) 

where A m is determined by Eq. (41) with the substitution M m. \5E\ obeys Eq. (50), 
L m is given by the angular momentum of the incident proton, and Lm is obtained by the 
alpha particle from the nucleus. We do not distinguish in a m = z$e 2 between zq and (z — 2) 
as the correction is of the same order as the small interaction. 

Since the energy transfer \SE\ is determined by small \ f—f\ ~ x , it strongly depends on 
angular momenta Lm and L m which set a smallest distance between particles. Lm and L m 
should be chosen to get the energy transfer \8E\ obeyed Eq. (50). We discuss this below. 
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E. Phase connection 



The action A relates to the phase connection 



^ 2 ~ A exp - 



A 



(58) 



of two states, (i) and (2), shown in Fig. 10, where (2) is the physical final state of the 
reaction. The state (i) consists of the nucleus gj^Th with the alpha particle and proton 
close to it. (1) is the physical initial state with the incident proton. The phase connection 
between the states (i) and (1) reads 



ipi ~ ipi exp 



:A m {e, L m ) 



(59) 



The dominant contribution to the imaginary phases comes from motion in the Coulomb 
field. By means of Eqs. (58) and (59) one can obtain the probability W of the reaction 



235 
92 



in the form 



where 



W 



w 2 



IfTh + a + p 



~ exp —A 



A = A M (E - \SE\,L M ) + A m {e + \6E\,L m ) - A m (e, L m ) 



(60) 



(61) 



(62) 



This procedure of phase connection is same as in III. 

One can use for the action (62) the approximation (41), which is used in the derivation 
of Eq. (50). In this approximation A does not depend on L m due to the cancellation in the 
last two terms of Eq. (62). This means, that L m can vary to adjust the energy transfer \5E\ 
to one given by Eq. (50) with no impact on A. In this situation the minimum A is reached 
at L M = 0. Finally, the action in Eq. (62) takes the form 



A = 



2M 



h VE-\5E\ 



-4 



X 



h 2 /2Ma M 



+ 



'2m 



h 



y/e + \5E\ 



~ -r, (63) 



where \SE\ satisfies the relation (50). 

The physical trajectories of particles in the nuclear reaction (60) are shown in Fig. 11. 
The probability of the channel (a) in Fig. 11 is almost 100% and the probability of the 
cannel (6), which is the reaction (60), is of the order of exp(— A) « 1. In Fig. 12 the 
classical positions of particles are shown at the moment of time when the incident proton 
reaches its minimum distance to the nucleus. 

The angular momentum L m of the proton should be exactly of the value to provide 
the energy transfer (50). Otherwise, the proton and the alpha particle do not meet at the 
nucleus position at the coincident "times" tm = r m , their interaction would be small and the 
action would be not a result of a cooperative motion of two particles but simply conventional 
A M (E, L M )- 
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F. Required angular momentum of the incident proton 



Let us calculate the energy transfer \8E\ at zero angular momentum of the incident 
proton L m = and then find the proton energy e which corresponds to this process. Since 
L m = 0, only re-components are involved which are determined by the classical dynamical 
equations 

= + (r x - R x f ; m ^ ~ {r x -R x y l ° } 

close to r (t ~ r M ~ r m ) the solutions have the form 



R x (ir) rMr) 



2/3 



Rs T s \ T m 

where R s and r s are some constants. The energy 5E, gained by the alpha particle, 



(65) 



P dr dR x 

SE = a >">l Tr^I* (66) 

diverges close to r and should be cut off by the condition R s (l — r/r m ) 2 / 3 > x . This gives 

i^i ^{wr'T (67) 

The ratio R s /r s , as one can see after a little algebra, satisfies the relation 



Ma m ( \L x ' 



ma M \r s 



I _ ^ s \ 2 _|_ aint 



l _R 1 V _a^fR£\ 2 (6g) 

r s J a M \r s J 



Substituting parameters for the reaction (60) M/m = 4, auloL m = 2, and a.i nt /a.M = 1/90, 
one can obtain R s /r s ~ 0.715 and the energy transfer \5E\ ~ 1.89 MeV. We use the estimate 
(37) for the nuclear size. With this \SE\ and the alpha particle energy E = 4.678 MeV, the 
relation (50) gives an unphysical (negative) value of e. This means, that an incident proton 
with zero angular momentum transfers too big energy and the real process requires a finite 
angular momentum L m in order to effectively increase the minimum distance xq in Eq. (67) 
to reduce the energy transfer. The minimum proton-nucleus distance increases when the 
proton centrifugal energy becomes of the order of the Coulomb one at the nucleus radius. 
This corresponds to the estimate (40) for proton 



L 2 a 



rn 



2mx1 Xq 

and a typical angular momentum of the incident proton should be L m ~ 10h. 



(69) 
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G. Euclidean resonance 



If to insert the energy transfer \5E\ from Eq. (50) into Eq. (63), one can obtain 




At e — e. 



where 




1/3 



e. 



max 



E 



~ 1.85 MeV 



(71) 



(E = 4.678 MeV) the energy transfer \8E\ = 0, the parameter L m /2mx a m has its maximum 
value, and the action (70) matches the conventional one (41) resulting in the tunneling 
probability 



The result (72) reasonably describes experimental data if to multiply (72) by the nuclear 
attempt frequency 10 21 s _1 . 

At e < Emax the energy transfer \5E\ becomes finite, L m decreases, and the action 
(70) reduces compared to Am(E,0). Upon reduction of e, the action (70) turns to zero at 
the certain proton energy £ R , which relates to Euclidean resonance. For the reaction (60) 
e R = 0.25 MeV and the accompanied energy transfer is \8E\ = 1.15 MeV. In other words, 
when in the reaction (60) the energy of the incident proton is e = 0.25 MeV, it converts 
into the proton with the energy e + \5E\ = 1.40 MeV and the energy of the emitted alpha 
particle (instead of E = 4.678 MeV) becomes E - \5E\ = 3.53 MeV. 

The cross-section of the reaction (60) is not exponentially small at £ — Er, it has a sharp 
peak at this proton energy, and is determined by the angular momentum L m which provides 
the energy transfer \5E\ = 1.15 MeV. A rough estimate of this angular momentum (69) 
results in the impact parameter h ~ 8xo, according to the classical relation L m = h\ / 2m£. 
A geometrical estimate of the cross-section at £ — E max is reduced to the ring area o ~ l-nhbh 
of the width 5h near the circle of the radius h, as shown in Fig. 12. Since h is some optimum 
value, related to an extreme action, Sh can be estimated as Sh/h ~ \j\J A M . By means of 
the relation (72), the geometrical estimate of the cross-section produces 15 nuclear impact 
areas a ~ 15(7TXq). 

At £ — £ m ax the exit point of alpha particle is 7.3rr and one of the proton is 10a;o. The 
incident proton is stopped by the nuclear Coulomb field at 54x , as shown in Fig. 11. The 
interaction of alpha particle with a moving proton is analogous to its interaction with some 
non-stationary field, which results in Euclidean resonance. As in Fig. 6, the exponentially 
small part of alpha particle wave function merges at some moment of time the growing 
part of the wave function. Trajectories in imaginary time is only a convenient language to 
describe the effect. In real (physical) time alpha particle does not approach the nucleus and 
interacts with it solely via the Coulomb field. For example, this practically excludes spin 
interaction between them. 

The above calculations hold for spherical nuclei. Alpha emitters may be not spherical, 
but real parameters of nonsphericality unlikely essentially violate the estimates. 



W ~ exp [-A M (4.678 MeV, 0)] 



-80.75 



~ 10 



-35 



(72) 
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VII. CONCLUSION 



When a proton approaches a nucleus, which is an alpha emitter, it creates a nonstationary 
Coulomb interaction with the tunneling alpha particle. At the certain proton energy there 
are conditions for Euclidean resonance and the Coulomb barrier becomes transparent for 
the passage of the alpha particle. Normally, H 5 U emits alpha particle of the energy 4.678 
MeV. When the energy of the incident proton is close to its resonant value 0.25 MeV, it 
reflects with the energy 1.40 MeV and simultaneously the alpha particle is emitted with the 
energy 3.53 MeV. Beams of 0.2 MeV-scale protons are "cheap" since they can be technically 
created in a relatively easy way. For this reason, low energy protons can be used for practical 
applications, for example, in disactivation of the nuclear waste. 

The analytical calculation of Euclidean resonance on the semiclassical basis is given in 
Ref. [25]. In principal, one can solve numerically the initial Schrodinger equation for decay 
of a metastable state in presence of nonstationary field using the proper algorithm and ac- 
counting the boundary conditions [30]. However, there is a serious problem in such numerical 
calculations. As one can see from Fig. 6, the branch (3), related to Euclidean resonance, 
starts to form from the exponentially small branch (2) which is of the order of exp(— Awkb) 
at the well position. In order to get a numerical calculation accounted this effect, one should 
choose very short steps At in time (the number of steps is proportional to (At)" 1 ) and Ax 
in coordinate (the number of steps is proportional to (Ax)" 1 ). They have to satisfy the 
condition At ~ (Ax) 2 < exp(— Awkb), otherwise the effect would be missed. Suppose the 
numerical computation with 10 3 steps in time and 10 3 steps in coordinate requires one sec- 
ond. One can easily estimate the total computation time as exp(1.5AwKB) x 10" 10 days. 
Since for alpha decay A WK b — 80, the total computation time should be 10 42 days. For this 
reason, numerical computation should start not with the initial Schrodinger equation but 
with some semiclassical approach. 
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FIGURES 



V(x) 




FIG. 1. The path of tunneling is denoted by the dashed line. E is the energy of the metastable 
state, V is the barrier height, and a is the typical potential length. 
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V(x) 




FIG. 2. The particle can absorb a quantum and tunnel in a more transparent part of the barrier 
with the energy E + Kl (a). The process of the multiquanta absorption with the subsequent 
tunneling (b). 



22 




FIG. 3. The multiquanta emission with the subsequent tunneling in a less transparent part of 
the barrier. 
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FIG. 4. (1) is the initial particle state with the energy E and (2) is the final state of the energy 
E — \ 5E\ after tunneling. The state (i) has the same energy as (2). There are two ways to connect 
phases of the states (i) and (1): by quanta absorption (a) and by quanta emission (b). 
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V(x) 




FIG. 5. The electron emission from the metal (left) to the vacuum by applying the constant 
electric field which is weakly modulated in time. 
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FIG. 6. The branches of the wave function ip ~ exp(— ImS). In the absence of a nonstationary 
field there are the dominant branch (1) and the sub-dominant one (2) which merge at the point 
"exit" and convert into the outgoing wave, denoted by the dashed line. With a nonstationary 
field, at t = the branch (3) is generated from the circled region at the sub-dominant branch (2) 
where the semiclassical condition is violated at t = 0. The branch (3) is formed during a short 
(non-semiclassical) time and then it moves semiclassically, keeping its maximum at the classical 
trajectory point x c i(t). This maximum value of the branch (3) decreases slowly in time due to 
quantum effects of smearing. 
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potential 




FIG. 7. The potential barrier which is passed by the alpha particle. The nuclear Coulomb field 
is cut off at the nuclear radius xq. R e is the classical exit point. 



27 



FIG. 8. The classical trajectory of the proton moving in the nucleus Coulomb field. r e is the 
classical turning point. 




FIG. 9. The classical trajectories in imaginary time of the proton r and of the alpha particle 
R. The both particles occur at the nucleus position xq at the same "moment" of imaginary time. 
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FIG. 10. (1) is the initial state of the reaction which includes the uranium nucleus and the 
incident proton. (2) is the final state, including the outgoing alpha particle and the proton. The 
state (i) consists of the thorium nucleus with the alpha particle and the proton close to it, but 
away of nuclear forces. The state (i) serves to connect the phases of the states (1) and (2). 




29 



FIG. 11. The big circle is the nucleus. The incident proton (the smallest circle) moves along 
the thick curve and may reflect elastically with no stimulation of alpha decay (a). The incident 
proton can initiate the alpha decay and then it appears, together with the alpha particle (the small 
circle), in the channel (b). The exit points at the x-axis are given for the energy of the incident 
proton 0.25 Mev. 




FIG. 12. h is the impact parameter of the incident proton. The cross-section of the reaction is 
determined by the area between two dashed circles, separated by the distance Sh. 
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